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ABSTRACT 

We report high spatial resolution 11.2 and 18.1 /ini imaging of V838 Monocerotis obtained with Gemini 
Observatory's Michelle instrument in 2007 March. Strong emission is observed from the unresolved stellar 
core of V838 Mon in our Gemini imagery, and is confirmed by Spitzer MIPS 24 /xm imaging obtained 
in 2007 April. The 2007 flux density of the unresolved mid-infrared emission component is ~2 times 
brighter than that observed in 2004. No clear change in the net amount of 24/im extended emission is 
observed between the 2004 and 2007 epoch Spitzer imagery. We interpret these data as evidence that 
V838 Mon has experienced a new circumstellar dust creation event. We suggest that this newly created 
dust has condensed from the expanding ejecta produced from V838 Mon's 2002 outburst events, and 
is most likely clumpy. We speculate that one (or more) of these clumps might have passed th rough 
the lin e-of- sight in late 2006, pro ducing the brief multi- wavelength photometric event reported by iBondl 
(|2006f ) and lMunari et all ()2007bf ). We detect no evidence of extended emission above a level of ~1 mJy 
at 11.2 /xm and ^7 mJy at 18.1 /im over radial distances of 1860 - 93000 AU (0'.'3-15'.'0) from the central 
source. Using the simple assumption that ejecta material expands at a constant velocity of 300-500 km 
s~^ , this gap of thermal emission suggests that no significant prior circumstellar dust production events 
have occurred within the past ~900-1500 years. 

Subject headings: circumstellar matter — stars: individual (V838 Monocerotis) — 



1. INTRODUCTION 

V838 Monocerotis experienced three photometric out- 
bursts in early 2002, beginning with an initial event 
on 2002 January 6.6 (jBrown et al.l l2002f ) and contin- 
uing with events i n 2002 February and 2002 March 
(jKimeswenger et al.l 2002; Munari et al. 2005). The star's 
spectrum cooled throug h the F, G, and K spectral types 
throughout February (|Munari et al.ll2007a| ), transitioned 
to a late M-type giant by 2002 April 16 (jRauch et al.l 
|2002| ). and had cooled to a L-type supergiant appear- 
ance by 2002 October (JEvans et al.l [2003il . The circum- 
stellar ejecta from the outburst events was not spherically 
symmetric, as revealed by the intrinsic polarization de- 
tected until 2002 Fe b ruary 13 jW isniewsM^JalJ l2003at 
lDesideraet"aIll2004[ ): IWisniewskij (|2007ll estimated that 
ejecta located interior to the region producing Ha and Ca 
II emission was flattened by a minimum of 10% from a 
spheroidal shell. A renewed intrinsic polarization compo- 
nent observed on 2002 October 22 might be attributable 
to a new source of asymmetrical scatterers in V838 Mon's 
circumstellar enve lope or to a change in the opacity of the 
existing envelope (jWisniewski et al.ll2003b[ ). 

V838 Mon produced a spectacular li ght echo fi rst de - 



source of this emitting dust could have both interstellar 
and circumstellar components. 

V838 Men resid es in a young open cluster 
(Afsar & Bond' '2007D at a distance of 6.2 ±1.2 kpc 
( Sparks et al.. .20081) . It has a known, albeit unresolved, 
B3V binary companion (jMunari et al.l |2005[ ) . The na- 
ture of V838 Mon prior to outburst is a subject of active 
debate (c.f. iMunari et all 120051: PTyfenda et all [2005'): a 
wide variety of mechanisms have been proposed to explain 
the eruptive even ts in V838 Mon, inc luding a shell ther- 
monuclear event (jMunari et al.ll2005[ ). a born-again star 
which experience d an accr etion event (jLawlorl l2005f ). a 
stellar merger (Tvlenda fc Soker 200a), and a star which 
engulfed several planets (JRetter et al.ll2006D . Determin- 
ing the amount and spatial distribution of dust within 
the V838 Mon line of sight, and quantifying how much 
is circumstellar versus interstellar in nature, would help 
constrain whether its outburst mechanism was a variant of 
a repeatable phenomenon, like a thermonuclear runaway 
or nova- like explosion, or more likely a singular event, like 
a stellar merger. 



tected on 2002 February 17 by Henden et all ( 2002 ) 



which evolved dramatically over time (jBond et al.l 



2003 ) 



2. OBSERVATIONS 

We obtained N (11.2^m) and Qa (IS.l^m) imaging of 
V838 Mon on 2007 March 21-22 using Gemini Oberva- 
tory's Michelle imager. Michelle is a 320 x 240 pixel ar- 
ray which has a 32'.'0 x 24'.'0 field of view (O'.'l per pixel). 
All observations were obtained using a standard chop-nod 
mode with the same instrument configuration: a detector 
position angle of 110° east of north, a chop angle of 130° 
east of north, and the maximum chop throw of 15'.'0. The 
weather and seeing conditions during our observing runs 
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iBaneriee et"all ()2006D imaged V838 Mon with the Spitzer 
Space Telescope at 24, 70, and 160/xm, and attributed the 
observed excess of thermal infrared (IR) emission from the 
unresolved central source as evidence of the formation of 
circumstellar dust produced by the 2002 outbursts. An 
IR excess was also observed f r om th e spatially resolved 
fight echo, and IBaneriee et all ()2006t ) concluded that the 



were stable, with FWHM values of individual chop-nod 
pairs deviating from the average values quoted in Table [1] 
by less than 0'.'02 on 2007 March 21, 0'/05 on 2007 March 
22, 0'.'04 on 2007 August 6, and O'.'Ol on 2007 August 9. 

The data, summarized in Table 1, were reduced with 
procedures similar to those prescribed in Gemini's midir 
package. Following combination of individual chop-nod 
pairs, known image artifacts such as vertical and horizon- 
tal "staircases" were characterized by sampling 15 pixel- 
wide regions along column and row borders, and removed 
via custom IDL routines. Cleaned images were regis- 
tered to a common position and examined for evidence of 
anomalous elongation along the chop direction, which is 
known to characterize a fraction of all Michelle data since 
early 2007. Our final combined images were flux calibrated 
using observations of Cohen standards (jCohen et al.lll999l ) 
HD 66141 and HD 60522 obtained immediately following 
our V838 Mon observations. Figure [1] provides a 15'.'0 x 
15'.'0 view of our reduced V838 Mon and PSF-star data, 
plotted on a linear scale in units of mJy. 

V838 Mon was also imaged w ith the Multiband Imaging 
Photometer for Spitzer (MIPS; iRieke et al. 2004) on 2007 
April 10 by program 30472 (PI K. Su). These data were 
reduced using version 16.1.0 of the Spitzer pipeline. Ear- 
lier epoch MIPS imaging from 2004 October 14 (24/xm) 
a nd 2005 April 2 (70/im) have been previously presented 



in lBaneriee et al.l (|2006) . 



3. ANALYSIS 
3.1. Photometry of the Unresolved Point Source 

We used simple aperture photometry to extract flux 
densities for the unresolved V838 Mon point source in 
our Gemini data. We explored the use of 5 circular aper- 
tures having radii ranging in size from 0'.'6-3'.'0 for our N 
data and from 0'.'7-3'.'0 for our Qq data, and found consis- 
tent results from each. The mean flux density values from 
these apertures, 29.54±0.14 Jy at 11.2/zm and 36.43±0.82 
Jy at 18.1/im, are cited in Table [T] The quoted errors 
are the standard deviations from our different apertures, 
and provide an estimate of the internal uncertainty in our 
measurements (~1%). Additional uncertainties in the ab- 
solute flux calibration, owing to the internal accuracy of 
the Cohen flux standards and changes in the atmospheric 
conditions between our observations, are expected to be of 
order several percent for bright sources observed in stable 
conditions (K. Volk 2008, personal communication). 

Spitzer's MIPS observed V838 Mon in 2007, --2 weeks 
after our Gemini data were obtained. The 24/zm im- 
agery was heavily saturated at the core. iBaneriee et al.l 
(|2006[ ) reported similar saturation effects in their 2004 
24/zm MIPS data. These authors fit a scaled PSF-star 
to match the brightness of V838 Mon's first airy ring, and 
replaced the saturated V838 Mon core by the core of this 
scaled PSF-star to recover a flux density measurement of 
the unresolved V838 Mon central source (Tabled]). To ex- 
tract similar measurements for the 2007 24^m MIPS data, 
we simply scaled and registered the 2004 and 2007 V838 
Mon data to match their unsaturated PSF wing regions. 
The optimal scaling factor yielded a 24/im flux density of 
29.69± }|f Jy (Tabled]). Note that the quoted error rep- 
resents the uncertainty we estimated for our image scaling 
(~1.5%) along with the 5% absolute flux uncertainty of 



the 20 04 Spitzer flux densities reported bv IBaneriee et all 
(|2006t ). Using simple aperture photometry, we also ex- 
tracted a flux density for the 2007 April Spitzer MIPS 
lOfim observation of 7.34±0.73 Jy (Table [1]). 

The 2007 March Gemini 18.1^m flux density is consis- 
tent with the 2007 April Spitzer 24^m flux density at the 
3ct level; these flux densities are a factor of two stronger 
than that observed in 2004 October observations. Simi- 
larly, the 2007 April Spitzer 70/im flux density is ~2 times 
stronger than that observed in 2005 April. The amount 
of excess IR emission associated with the unresolved cen- 
tral source of V838 Mon has clearly increased between 
2004/2005 and 2007. 

3.2. Extended Emission 

Gemini's Michelle provides a 10-fold improvement in 
spatial resolution compared to Spitzer's MIPS; hence, 
our data provide us the unique opportunity to search 
for mid-IR extended emission interior to that probed by 
IBaneriee et al.l ()2006D . Comparison of median azimuthally 
averaged radial profiles of our 11.2/xm and 18.1/Ltm im- 
agery of V838 Mon with a similar color PSF-star obser- 
vation (HD 52666; M2 III), obtained immediately before 
our V838 Mon data, reveal no convincing evidence of ex- 
tended emission nearby the stellar core. Measured median 
FWHM values of V838 Mon and our PSF-star are the same 
at 11.2 and 18.1 /im (Table]!]) to within the maximum ob- 
served deviation amongst chop-nod pairs (<0'.'02 at 11.2 
/im and <0'.'05 at 18.1 /im), which quantitatively supports 
this conclusion. 

To search for extended emission outside of the imme- 
diate vicinity of the stellar core, we used both contem- 
poraneously observed PSF-star observations (HD 52666; 
Table [1]) and longer integrations of a PSF-star obtained 
several months after our V838 Mon observations (HD 
168723; Table[l]) to model and subtract the PSF from our 
V838 Mon imagery. The deeper PSF-star observations 
exhibit marginally different FWHM values as compared 
to V838 Mon, likely attributable to minor differences in 
atmospheric conditions, but provide a lower average de- 
tector noise background which is better suited to iden- 
tify faint extended emission far from the V838 Mon stellar 
core. Figure [3] illustrates our V838 Mon imagery following 
subtraction of optimally scaled and registered PSF-star 
observations, plotted on a linear stretch in units of mJy, 
and smoothed by a second-order Gaussian function. We 
have masked the region inside a radial distance of 0'.'8 (8- 
pixels), which is dominated by PSF-subtraction residuals, 
along with a diagonal strip in our N imagery which ex- 
hibited a commonly observed detector artifact which arises 
when Michelle images bright sources. 

We detect no convincing evidence of extended emission 
above a level of ~1 mJy at 11.2 /im and above a level of 
~7 mJy at 18.1 /im. Applying additional binning or alter- 
nate smoothing functions to the data presented in Figure 
[3] yielded no evidence of extended emission above the ob- 
served residual detector noise background. 

As described in Section 13.11 we shifted and scaled the 
archival 2004 Spitzer MIPS 24/im data to extract point 
source photometry for the 2007 Spitzer 24/xm data. The 
scaling factor we applied to the 2004 data (1.97) led to 
an optimal cancellation of the PSF wings of the 2007 



data in the subtraction process, but produced a large re- 
gion of over-subtracted flux exterior to the stellar core, 
whose morphology correspo nded to th e extended emis- 
sion reported by iBaneriee et al.1 (|2006[ ). Using a circu- 
lar aperture of radius Sff.'O and masking the central PSF- 
core region, we determined the net 24/im flux density of 
this oversubtraction region was -0.88 Jy. This is exactly 
the flux density one would expect to see if the extended 
emission component did not vary between 2004 and 2007, 
;iven the reported 20 04 extended emission flux of 0.9 f Jy 
Baneriee et al.l |2006[ ) and our use of a scaling factor of 
1.97, i.e. 

24^7712007 " (1-97 X 24^7^2004) = 24/i777„et 

0.91Jy - (1.97 X 0.91Jy) = -O.SSJy 

While a comprehensive comparative analysis of the 2004 
and 2007 Spitzer MIPS data is outside of the scope of 
this Letter, our simple analysis suggests a) the amount 
of 24/zm and 70/im emission in the unresolved V838 Mon 
stellar core has increased by a factor of ^2 from 2004 to 
2007; and b) the net amount of 24/im extended emission 
appears unchanged from 2004 to 2007. 

4. DISCUSSION 

We detect no evidence of extended 18.1 /im emission 
above a level of ~7 mJy from radial distances of 1860 AU 
0"3) to 93000 AU (15^/0), assuming a distance of 6.2 kpc 
Sparks et al.l 120081 ). Making the simple assumption that 
V838 Mon's circumstellar ejecta expands at a constant ve- 
locity, this gap of thermal IR emission suggests no circum- 
stellar dust producing events h ave occurred within the p ast 
900 {ifvejecta ~300 km s~^: IWisniewski et al.ll2003aD t o 



~1500 years (if Vejecta ~500 km s~^: lMunari et al.ll200l 
Our analysis of 2007 imagery of V838 Mon from both Gem- 
ini's Michelle and Spitzer's MIPS indicates the presence of 
significantly enhanced 10-70/im emission in the unresolved 
central stellar core in 2007 as compared to 2004. We in- 
terpret our 2007 data as evidence that a new circumstellar 
dust production event has occurred since 2004. 

Spectroscopic evidence of V838 Mon's binary compan- 
ion disappeared and Ha returned to emission in late 2006, 
corresp onding with a temporary m ulti-colo r photometric 
fading jBondl 120061: [M^nari et al.|[2007bD . iMunari et all 
(|2007bl ) attributed these events to either a) an eclipse be- 
tween the primary and secondary; or b) a dust cloud char- 
acterized by Eb-v = 0.55 and Ry = 3.1, gravitationally 
bound to the binary system, which eclipsed the B3V sec- 
ondary. We suggest alternate interpretations are possible 
if these photometric and spectroscopic events are corre- 
lated to the new mid-IR behavior we have found in data 
from early 2007. A gravitationally bound dust cloud would 
radiate a constant amount of thermal IR emission through- 
out its orbit except during the brief time that it is occulted 
by the secondary; thus, unless it was only recently formed, 
such a cloud could not be singularly responsible for pro- 
ducing the observed enhanced mid-IR emission. A simple 



eclipse between the primary and secondary also would not 
le ad to e nhanc ed thermal IR emission. 

iBondl (|2006r ) suggested that the photometric and spec- 
troscopic events of late 2006 were indicative of the 2002 
ejecta interacting with and overtaking the B3V binary. 
More recent observations indicate that V838 Mon is expe- 
riencing a prolonged period of reduced photometric bright- 
ness (H. Bond 2008, personal communication), further 
supporting this interpretation. While the temporal res- 
olution of our data do not allow us to quantitatively de- 
termine whether this proposed encounter is related to the 
new circumstellar dust creation event we report, we sug- 
gest that such a correlation is plausible. 

We suggest that the expanding ejecta from the 2002 
outbursts has condensed to form new circumstellar dust, 
producing the enhanced IR excess observed inside a ra- 
dial distance of 1860 AU. The dust condensation radius 
of ^3.5 AU for the B3V compani on, which assumes Tcond 
-lOOOK, Rcond ~ ^^(t^)^ (|lreland et al.ll2005[) . and 
the estimated primary-secon dary separation o f ^^34 AU 
derived from interferometry (jLanet et al.ll2005l ). indicates 
that dust can survive within the radial distance constraint 
of <1860 AU implied by our observa tions. This type o f 
event is not altogether unexpected, as I Lynch et al.l (|2004() 
previously reported the detection of an expanding region of 
molecular gas which they noted was likely dust precursor 
material. The enhanced flux at 18.1, 24, and 70 /xm we 
observed in 2007 suggests that the equally strong (2007 
epoch) 11.2 //m flux might not represent an additional 
enhancement of gas-phase molec ular emission above that 
suggested by iLynch et al.l (|2004l ). but rather represents a 
warm thermal dust emission component. 

Early post-outburst polarimetric ob servations 
(|Wisniewski et al.l l2003al : iDesidera etllj 12004 ^ indicated 
this ejecta was asymmetric; hence, we suggest that it is 
likely that dust which condenses from this ejecta will be 
spatially non-uniform and/or clumpy. We speculate that 
it is possible that one or more of these clumps eclipsed the 
line of sight in late 2006, pro ducing the attenuation event 
reported by iBondl (|2006.) and lMunari et all (l2007bD . Un- 
like the dust cloud eclipse event invoked by iMunari et al.l 
(|2007bD . we believe that any eclipsing dust blob is more 
likely to located in the expanding ejecta which is con- 
densing to form new dust. FoUowup interferometric ob- 
servations might reveal differences in the angular size of 
the system from that found in 2004 epoch observations 
by iLanet et al.l (|2005f ) if this ejecta has indeed recently 
reached the location of the B3V secondary star. 
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Table 1 
Summary of Observations and Flux Densities 



Target 


UT Date 


Observatory 


^central 


<5A 


FWHM 


Flux Density 


Comment 








(^m) 


(^m) 


(arcsec) 


(mjy) 




V838 Mon 


2004 October 14^ 
2005 April li 


Spitzer 
Spitzcr 


23.7 

71.42 


4.7 
19 


O'.'O 
18"0 


15.06±0.75^ 
3.82±0.382 


Baneriee et al. (2006) 


V838 Mon 


Baneriee et al. (20061 


HD 52666 


2007 March 21 


Gemini 


11.2 


2.4 


0'.'40 




PSF-star obs 


HD 66141 


2007 March 21 


Gemini 


11.2 


2.4 


0'.'43 




flux std 


V838 Mon 


2007 March 21 


Gemini 


11.2 


2.4 


0'.'40 


29.54±0.143 




V838 Mon 


2007 March 22 


Gemini 


18.1 


1.9 


0'.'61 


36.43±0.823 




HD 52666 


2007 March 22 


Gemini 


18.1 


1.9 


0'.'57 




PSF-star 


HD 60522 


2007 March 22 


Gemini 


18.1 


1.9 


0'.'59 




flux std 


V838 Mon 


2007 April 10 


Spitzer 


23.7 


4.7 


O'.'O 


29.67±;|52 


AOR-21475072 


V838 Mon 


2007 April 10 


Spitzer 


71.42 


19 


IS'.'O 


7.34±0.732 


AOR-21475072 


HD 168723 


2007 August 6 


Gemini 


18.1 


1.9 


0'.'44 




deep PSF-star obs 


HD 168723 


2007 August 9 


Gemini 


11.2 


2.4 


0'.'54 




deep PSF-star obs 



Note. — ^ 2004/2005 Spitzer flux density measurements were adopted from Table 2 of lBaneriee et all |)2006I V ■^ The uncertainties quoted for 
Spitzer flux density measurements, 5% for 2Afj,m and 10% for 70^m, represent the estimated absolute flux uncertainty. ^ Gemini flux density 
error estimates represent the internal uncertainties of the measurements; additional absolute flux scaling uncertainties are expected to be of 
order several percent for bright sources observed in stable conditions (K Volk 2008, personal communication). 




Fig. 1. — Reduced imagery spanning a 15''0 x 15'.'0 field of view and plotted on a linear scale in units of injy. The displayed panels include 
(a) N'-band V838 Mon; (b) N'-band HD 52666 (PSF-star); (c) N'-band HD 168723 (deep PSF-star); (d) Qa-band V838 Mon; (e) Qa-band 
HD 52666 (PSF-star); and (f) Qa-band HD 168723 (deep PSF-star) imagery. The diagonal dark lanes seen prominently in panels a-d are 
known detector artifacts which arise when Gemini's Michelle observes bright sources. 
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Fig. 2. — SED of the unresolved point source flux of V838 Mon, based on the flux densities listed in Table[T] The 2007 Gemini and Spitzer 
data, obtained ~2 weeks apart, are consistent to within Scr of the total error budget of each observation. The IR excess emission at ~20/.tm 
and 70/xm has clearly increased between 2004 and 2007. We interpret this as evidence that a circumstellar dust production event has recently 
occurred. 




Fig. 3. — PSF-subtracted V838 Mon imagery plotted on different linear scales in units of injy. The field of view, IS'.'O x IS'.'O, lies entirely 
inside of region of extended emission detected by Spitzer's MIPS in 2004 l IBaneriee et al.ll2006l ). The displayed panels include (a) N -band 
imagery using HD 52666 as the PSF-star; (b) N -band imagery using HD 168723 as the PSF-star; (c) Qa-band imagery using HD 52666 as the 
PSF-star; and (d) Qa-band imagery using HD 168723 as the PSF-star. Regions inside a radial distance of 0'.'8 (8-pixels) were contaminated 
by PSF subtraction residuals and are therefore masked in the displayed imagery. We detect no clear evidence of excess extended emission in 
these data. 



